A single cold atom as efficient stationary source of EPR-entangled light 
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The Stokes and anti-Stokes components of the spectrum of resonance fluorescence of a single 
trapped atom, which originate from the mechanical coupling between the scattered photons and 
the quantized motion of the atomic center of mass, exhibit quantum correlations which are of two- 
mode-squeezing type. We study and demonstrate the build-up of such correlations in a specific setup, 
which is experimentally accessible, and where the atom acts as efficient and continuous source of 
EPR-entangled, two-mode squeezed light. 
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I. INTRODUCTION 

The control of atom-photon interaction is object of 
intensive research for its potentialities in quantum net- 
working. In fact, several experimental realizations have 
accessed novel regimes of engineering atom-photon in- 
teractions and have opened promising perspectives for 
implementing controlled nonlinear dynamics with simple 
quantum optical systems. Fundamental steps in this di- 
rection have been, amongst others,the generation of en- 
tangled light in atomic ensembles [ILI4L1 atomic memory 



for quantum states of light |3l, 4, [5], |6| , 0], and entan- 
glement of remote ensembles [8|, [9|, |l0|. At the single 
atom level, entanglement between a single atom and its 
emitted photon [ll[ has been demonstrated in pjl, [l3| , 
while in cavity quantum electrodynamics generation of 
quantum light has been achieved, like lasing at the single 
atom level ^A, [lfij], controlled single-photon generation 
[3 E3, EH [l9 1) as well as quantum state and entangle- 
ment engineering in the microwave regime [20| . 

Quantum networking with single trapped atoms or ions 
shows several advantages, due to the high deg ree of con- 
trol one can achieve on these systems [ill [22| . Con- 
trol can be gained on the internal as well as on the 
external degrees of freedom, which can both be inter- 
faced with light by exchange of angular and linear mo- 
mentum. In particular, by coupling the atomic external 
degrees of freedom with photons via the mechanical ef- 
fect of light, atom-photon interfaces for continuous vari- 
ables can be imp lemented even at the level of a single 



atom 



applied in [25 



This concept has been specifically 
where the realization of a pulsed op- 
tical parametric amplifier based on a single cold trapped 
atom inside a high-finesse optical cavity was proposed, 
and it was shown theoretically that this system allows 
for the controlled, quantum-coherent generation of en- 
tangled light pulses by exploiting the mechanical effects 
of atom-photon interaction. 



In this manuscript we investigate the quantum corre- 
lations between the Stokes and anti-Stokes sidebands of 
the resonance fluorescence of a trapped atom, i.e. be- 
tween the spectral components which are due to the cou- 
pling of the electro mag netic field to the atom's oscilla- 
tory motion [27l. [28l. [29L f30l| . The spectrum is studied for 
an atom tightly confined inside a resonator and continu- 
ously driven by a laser, in the setup sketched in Fig. [1] 
This setup has been considered in [2^, [26[ for the case of 
pulsed excitation, where scattering could be considered 
coherent. In the present work, the atom is continuously 
driven and hence both coherent and incoherent scatter- 
ing processes determine the dynamics of the system. We 
find that in a suitable parameter regime the Stokes and 
anti-Stokes spectral components of the resonance fluo- 
rescence are two-mode squeezed, that is, their amplitude 
and phase quadratures are quantum correlated. In fact, 
the variance of the difference of the amplitude quadra- 
tures of the two sideband modes, as well as the variance 
of the sum of their phase quadratures, are squeezed be- 
low the shot noise limit, hence reproducing the salient 
properties of the entangled, simultaneous eigenstate of 
relative distance and total momentum of two particles, 
as considered in the original EPR paradox [3l|, |32 |. In 
our model, entanglement between the modes originates 
from the mechanical coupling of the electromagnetic field 
with the quantum motion of the atom, and it is endorsed 
by a specific setup, which achieves resonant emission of 
the Stokes and anti-Stokes photons. In this regime, the 
single atom acts as an efficient continuous source of EPR- 
entangled, two- mode squeezed light. 

Conventionally, two-mode squeezed states emerge from 
the nonlinear optical interaction of a laser with a crys- 
tal, i.e. from parametric amplification or oscillation. As 
such, the phenomenon is the result of many-atom dy- 
namics, often described by a simple nonlinear polariza- 
tion model. In the single-atom case novel features ap- 
pear which are due to the coherent microscopic dynam- 
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ics. Our study allows us to identify the dependence of 
these features on the external parameters, thereby giv- 
ing us insight into how macroscopic properties arise from 
microscopic dynamics in this particular non-linear pro- 
cess. Moreover, we find peculiar spectral characteris- 
tics of the squeezing which are unique to this system, 
and which we trace back to the interplay of the var- 
ious time scales of the dynamics. In a more general 
context, our study is connected to previous work on the 
quantum features of the spectrum of resonance fluores- 
cence [S3, US US H3, [H, Ha H3, [H, |H, S3] , and to recent 
experimental and theoretical studies on quantum correla- 
tions in the light scattered by atoms 0, 0, S 1, SH S3, El , 
by semiconductor microcavities [HJ , and by macroscopic 
mirrors [HIHl]. 

This article is organized as follows. In Sec.[TT]the basic 
coherent dynamics, giving rise to quantum correlations 
between the Stokes and anti-Stokes components of the 
spectrum of resonance fluorescence, are briefly reviewed, 
and the important time scales are introduced. In Sec. IIIII 
the theoretical model is described in detail and the rele- 
vant scattering processes in the system arc identified and 
discussed. In Sec. |TV] the spectrum of squeezing is eval- 
uated using Quantum Langevin Equations; for a quick 
overview of the main results without the full theoretical 
elaboration, the reader may first skip this part and jump 
to Sec. IIV CI where the squeezing characteristics are cal- 
culated for a specific, experimentally achievable physical 
system. Finally, Sec. |V] presents the conclusions and an 
outlook. 
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FIG. 1: Layout of the system. A single atom is confined by 
an external potential inside an optical cavity and is driven 
by a laser. The cavity is resonant with the motion-induced 
Stokes and anti-Stokes components of the resonance fluores- 
cence. Correlations between these spectral component are 
measured in the cavity output. The orientation of the con- 
sidered vibrational mode has non-zero projection onto the 
laser direction. A possible geometry to implement the sys- 
tem would be an F = to F' — 1 atomic transition with the 
quantization axis B along the cavity axis, and B, laser wave 
vector, and laser polarisation mutually orthogonal, and a mo- 
tional mode parallel to the laser direction. More details can 
be found in Refs. [25l . [2(| where pulsed coherent excitation 
was considered. In the present paper we deal with continuous 
laser excitation. 



that |x2 1 > |xi|) with an angular frequency 



O = V\X2\ 2 - \xi\' 



(2) 



II. EPR-ENTANGLEMENT OF LIGHT AT THE 
CAVITY OUTPUT 

In this section we briefly review the coherent dynam- 
ics, described previously in Refs. [ID, which lead to 
two-mode squeezing between the Stokes and anti-Stokes 
modes in the light scattered by a trapped, laser-driven 
atom. We thus first ignore incoherent processes and fo- 
cus on the pulsed dynamics which can be obtained in 
a suitable parameter regime with a setup like the one 
shown in Fig. [TJ 



The time-evolution of the operators, in the Heisenberg 
representation, is given by [46[ 
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The trapped atom is coupled to an optical cavity of which 
two modes are resonant with the Stokes and anti-Stokes 
sidebands, respectively. For short times the laser-induced 
resonant interaction between the center-of-mass oscilla- 
tion, denoted by annihilation and creation operators b 
and w, and the two cavity modes, represented by oper- 
ators dj and aj (j = 1,2), is described by the effective 
Hamiltonian in the interaction picture 



Wcff = i?ixi a I^ + iS-X2fl2^ + H.c 



(1) 



where the scalars Xj indicate the strength of the coupling. 
This Hamiltonian generates periodic dynamics, provided 



In general these solutions describe tripartite entangle- 
ment among cavity modes and center-of-mass oscilla- 
tor [11 . An interesting situation is found after half a 
period, for T v — tt/Q. At this time (modulus 27r) the 
center-of-mass oscillator is uncorrelated with the cavity 
modes, which exhibit EPR-type entanglement [25l. |26|. 

Clearly, this description is approximate, and valid only 
when incoherent processes can be neglected. In the 
present work we consider the situation in which the atom 
is continuously driven by the laser field, such that quan- 
tum noise and dissipative processes affect the dynam- 
ics relevantly. We show that steady state entanglement, 
i.e. quantum-correlated spectral fluctuations in the two- 
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mode cavity output field, is found also under these con- 
ditions. The details of this entanglement will depend on 
the comparison between the time scale set by the coher- 
ent dynamics, _1 , and the time scales of the dissipative 
processes, n^ 1 for loss of photons from the cavity, and 
7 _1 for spontaneous scattering from the atom. In partic- 
ular, we will show that the squeezing spectrum shows 
distinct, qualitatively different features in the regimes 
< k, = k, and > k. The reader is referred 
to Sec. IIV CI where the spectra for different parameter 
regimes are reported. 



III. SCATTERING PROCESSES 

The purpose of this section is to discuss the coherent 
and incoherent scattering processes determining the dy- 
namics of the system. We will present these processes 
using physical pictures derived from the scattering ma- 
trix under moderate simplifications, in order to illustrate 
the more rigorous derivations presented in the subsequent 
section. We first introduce the model, and then identify 
the scattering processes and determine the corresponding 
rates. 



A. Model 

We consider an atom of mass M inside an optical res- 
onator and driven by a laser. The atomic motion is con- 
fined by an external potential, which we assume suffi- 
ciently steep in the radial direction so that the motion 
in this plane can be considered frozen out. We denote 
by x the axis of the remaining one-dimensional atomic 
center-of-mass motion. Moreover, we assume that only 
the atomic dipole transition between ground state \g) 
and excited state |e) couples relevantly to the fields, such 
that we can restrict the electronic dynamics to these two 
states. The atomic dipole is laser-driven, and it cou- 
ples to two modes (j — 1,2) of the resonator, as well 
as to the external modes of the electromagnetic field. 
The cavity modes couple also to the external modes of 
the electromagnetic field through the imperfect mirrors 
of the resonator. The total dynamics is governed by the 
Hamiltonian 

H = H + W, 

where Ho is the self-energy of the system of atom and 
fields, and W describes their mutual interaction, as well 
as the coupling between the cavity modes and the exter- 
nal modes through the finite transmission at the cavity 
mirrors. We now introduce each term in detail, and dis- 
cuss the dynamics in the reference frame of the laser at 
the angular frequency u>l- We decompose Hq according 
to 

H = H a + H c + H emf . (6) 



Here, H a is the Hamiltonian for the relevant atomic de- 
grees of freedom, 

H a = -hA\e){e\+H mec , (7) 

where A = ujl — is the detuning of the laser from the 
dipole transition at the angular frequency too, and 

H aBB = nvUtb + ^ (8) 

describes the harmonic motion of the atomic center of 
mass at angular frequency v, as determined by an ex- 
ternal potential, where b, are the annihilation and cre- 
ation operators, respectively, of a quantum of vibrational 
energy Tiu. In particular, the atomic position is given by 
x = y/h/2Mv(b + We denote by \n) the eigenstates 
of -ffmec at energy fiv{n + 1/2). The Hamiltonian for 
the cavity modes, which couple appreciably to the dipole 
transition, is 

where 8j — uj^—uij are the detunings of the laser from the 
frequencies u>j of two optical modes, and Oj,aj are the 
respective annihilation and creation operators of a quan- 
tum of energy Tiujj , i.e. a photon in mode j. We denote by 
\n\, n%) the eigenstates of H c at energy —hSini — fibin-i-, 
and consider the situation in which the mode frequencies 
fulfill the relation 

oj 2 - wi = 2v' . (10) 

where 

i/ = v + 6v (11) 

and 8v takes into account radiative shifts, such that cav- 
ity modes 1 and 2 can be simultaneously resonant with 
the Stokes and the anti-Stokes transitions. This contri- 
bution will be discussed in Sec. lIIIB~4l and determined in 
Sec. ELS 

Finally, the modes of the electromagnetic field external 
to the cavity possess the free Hamiltonian 

ffcmf = -?iX^ k i r kj rk j ~ n ^Z 5 ^ r ij^ > 

kj k s 

where r\, r| are annihilation and creation operators, re- 
spectively, of a photon at angular frequency uj\ = uj^—dx, 
wavevector k>, and polarization e^. Here, the subscripts 
A = k s and A = kj indicate the modes of the field which 
couple to the dipole and to the cavity modes (through 
the mirrors), respectively. The interaction term 

W = H aL + H ac + M/ ks + W kj (12) 

describes the couplings among atom and fields, decom- 
posed into four terms which correspond to the coupling 
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between atom and laser (H a L), atom and cavity modes 
(H ac ), atom and modes of the external electromagnetic 
field (Wk s ), and cavity modes and external electromag- 
netic field (Wkj)- We discuss these terms in the Lamb- 
Dicke regime, when the atomic motion is well localized 
over the wavelengths of the fields, such that the Lamb- 
Dicke parameter rj — ^/Tik 2 /2Mv is small, 77 1 . At 
lowest order in », the coupling between laser and dipole 
has the form [47| 



H aL = nn^ 



l-+cos 2 i 



-(26*6+1)) (13) 



+ i7 7 cos0 L (& t + 6) + 0(t7 2 ) 



H.c. 



with a = \g){e\ the dipole lowering operator and <7+ its 
adjoint, £1 the Rabi frequency, and Qj, the angle between 
the direction of propagation of the laser and the motional 
axis x. In what follows we denote the moduli of all rele- 
vant wave vectors by A:, as their differences are negligible. 
The coupling between the dipole and the cavity modes is 
represented by 

,2 



H„ 



i=i,2 



1 



2 



3(26+6+1) 



T/cos0 c tan0j(6 t + 6) + H.c. + 0(7y 2 ) 



(14) 



where g« is the coupling strength of the dipole to mode j, 
and the cavity axis forms an angle 9 C with the axis x of 
the motion. The angle <j>j takes into account the position 
of the trap center inside the standing wave of the cavity. 
Finally, the terms 



k s cr'r ks 



2 



,(26+6+1; 



+ ir? cos 6^(6 + 6+) + 0(?7 2 ) +H.c), 
W kj = n 9k } {a]r kj + H.c.) 

describe the coupling of atom and cavity to the modes of 
the external e.m. -field. Here, Wk s is the coupling of the 
dipole, at Rabi frequencies g^ a , with the external modes, 
whose wave vectors form angles #k B with the motional 
axis. This coupling gives rise to the finite linewidth 7 of 
the excited state, 7 = 27rpk a (^o)l5k s (^o)| 2 , with pk s (wo) 
density of states of the e.m. -field coupling to the atomic 
dipole at angular frequency ujq. The term describes 
the coupling of the cavity modes with the external modes 
at strength . This coupling gives rise to the linewidth 
of the cavity modes Kj — 7r|<?kj | 2 Pkj (wj), with p^Juij) 
density of states of the e.m.-field coupling to the cavity 
modes at angular frequency uij. 



B. Basic scattering processes 



limit all terms of W are weak perturbations to the dy- 
namics. We assume that the system is in the initial state 



= |5,n;0i,0 2 ;0kj;0 ks ), 



(15) 



with energy Ei — Hun, where the atom is in the ground 
state \g), the center-of-mass oscillator is in the number 
state |n) , and the cavity modes and the external e.m.-field 
are in the vacuum state, |0i, O2; Ok-; 0k s ). The scattering 
matrix elements between the initial state and all possible 
final states \ipf), with energy Ef, have the form 



S if = % - 2mS(E f - Et)Tif 



(16) 



where 5,/ is the Kronecker-delta, S(Ef — Ei) is a delta- 
function giving energy conservation between initial and 
final states, and Tit is the transition matrix to be evalu- 
ated in lowest order in perturbation theory, 



E.; — H 



with 



H cS = -U (A + i~j |e)(e|+7w6 t 6-?i ^ ( S j + ' lK j) a ) a 3 

(17) 

We now consider all possible scattering transitions to res- 
onant states, i.e. to final states \ipf) at energy Ef = Ei. 



1. Scattering of laser photons into the external e.m.-field 

We consider the scattering of a laser photon 
into the external e.m.-field by spontaneous emission, 
hence coupling of \tpi) to the final states |^k„) = 
|g, n'; Oi, O2; Okj; lk s )- This process is sketched in 
Fig. HJa). Here, the coupling with the cavity mode is 
neglected, as the cavity is far-detuned from the dipole, 
and the rate of this process can be approximated by the 
scattering rate of the atom in free space, 
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(18) 



where 



" A + i 7 /2' 
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, (20) 



t^ = m n{ cos9l /o + -^-» . (21) 

- 1 " 1 A + v + ij/2 A + i 7 /2 ' 



We consider the limit in which the atom is far-detuned 
from cavity modes and laser, |A| 3> 7, 8j,gj,il. In this 



The process described by amplitude (119)) does not af- 
fect the dynamics of the cavity modes nor that of the 
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FIG. 2: Basic scattering processes, (a): A laser photon is absorbed and emitted by the atom, without coupling to the cavity 
mode, (b) A laser photon is scattered into the cavity mode and then rescattered by the atom into the external modes of the 
electromagnetic field, (c) A laser photon is scattered by the atom into the cavity mode, and then it is transmitted by the cavity 
mirror into the modes of the external electromagnetic field. 



center-of-mass motion. In contrast, the amplitudes ([20]) 
and (|2ip are coherent superpositions of scattering pro- 
cesses involving, respectively, the mechanical effect of 
the laser and of the emitted photon on the atomic mo- 
tion (28l . [3(3] , thereby affecting the coherence of the mo- 
tional state. Their rate is 76 » rj 2 (cos 2 9l + a)jQ 2 /A 2 , 
where a describes the angular dispersion of the spon- 
taneously emitted photons, determined by the quantum 
numbers of the atomic transition 14811 . 



2. Scattering of laser photons into the cavity modes 

Next we discuss the processes in which a laser pho- 
ton is scattered into one of the cavity modes, thereby 
coupling the initial state to the states \tpi) — 

\g,n'; li, 2 ; kj ; k J or \ip 2 ) = |ff,"';0i, l 2 ;0 kj ;0 ks ). As 
these states are not stable, but resonantly coupled to the 
continuum of states \g, n'; 0i, 2 ; l kj ; ks ) by cavity decay, 
the correct final states of these scattering processes de- 
scribe the processes sketched in Fig. [5Jc) and have the 
form 



1 



Qj 



£ ki -H 



w kj IV. 



(22) 



where Qj projects onto the subspace orthogonal to \ipj), 
and Z kj ensures the normalization of the state. Further- 
more, iT kj gives the occupation probability of state \ipj), 
since Z kj = |(^| V> kj ) ? ■ 

The coupling rate between state | an d states l^kj) 
takes the form 
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(24) 



t^ = riQgUos, 



= rjflg* cos < 
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A - v + i 7 /2 
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(2,5) 
J2.6) 



Like in Eqs. (|20|) and (j2lj) . we recognize on the RHS of 
Eqs. (|25|) and (|26|) the coherent addition of two scatter- 
ing amplitudes, here representing the mechanical effects 
of the laser and of the cavity, respectively [5fj, |5l| • These 
processes are at the basis of the coherent coupling be- 
tween the atomic motion and the cavity modes described 
by Hamiltonian |T]), where Xi — ~ i*i a + and \2 = — i*l a -- 
We are interested in the regime where energy can be 
stored in the cavity modes through this coupling, which 
requires IX2I ^ lb as a necessary condition. In this 
situation, it is visible from the equations that in the limit 
Kj <C v, by choosing 5\ = v and 82 = —v one can achieve 
the optimum enhancement of the scattering of a laser 
photon into mode 1 accompanied by the excitation of 
the motion by one vibrational quantum, and of the scat- 
tering of a laser photon into mode 2 accompanied by the 
de-excitation of the motion by one vibrational quantum. 
Note that these scattering terms, Vf±^ nave an incoher- 
ent component which scales with 7 /|A ± v\. Therefore, 
in general coherent dynamics can only be achieved when 
7 -C I A I, on a time scale such that incoherent terms are 
negligible. Moreover, the condition 7 -C v is also re- 
quired in order to create quantum correlations between 
the two cavity modes, since the difference between the 
two coupling strengths \i and xi determines the typi- 
cal time scale on which entanglement is established, see 
Sec. |n]and 



3. Scattering of cavity photons into the external e.m. -field 

Assuming that photons have been coherently scattered 
into the cavity modes, they can be re-absorbed by the 
atom and emitted spontaneously into the external e.m.- 
field, as sketched in Fig. (Hb). In order to focus on the 
evaluation of the corresponding element of the scattering 
matrix, we consider the regime of very small cavity loss 
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rate, i.e. we assume stable cavity modes and ignore, for 
the clarity of the picture, cavity decay. Be the initial 
state 



\1>i 



|5,n;mi,m 2 ;0 kj ;0 ks 



(27) 



at energy E i m = Tivn — hmiSi — Tim,2b~2, with the atom 
in \g), the center-of-mass oscillator in the number state 
\n), the cavity modes in the Fock states \m\) and |m 2 ), 
and the external e.m.-field in the vacuum state, |0 kj ; ks ). 
This state is coupled to the states 

l^/,mi> = |g,™; m i - l,TO 2 ;O kj ;lk s ) (28) 
\i>f,m' a ) = \g,n;mi,m,2 - l;0 kj ; l ks ) (29) 

by absorption of a cavity photon and spontaneous emis- 
sion. We evaluate the corresponding rate under the as- 
sumption, that tan</>j = 0, i.e., there are no mechanical 
effects of the resonator on the atom at first order in r), 
and find an effective loss rate of the cavity modes 
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A - Sj + i7/2 



(30) 



It should be noted that these processes arise from 
atomic scattering of a laser photon into the cavity modes, 
which is then rescattered by atomic emission into the ex- 
ternal modes of the e.m.-field. Hence, these processes 
can interfere with atomic scattering of a laser photon, 
in the limit discussed in Sec. IIII B II in which the cou- 
pling to the cavity plays no role. In these calculations we 
have not considered the coherent addition of these two 
noise effects, but we will consider phase relations and 
possible interference in these noise sources when study- 
ing the dynamics with the quantum Langevin equations 
in Sec. ITVAl 

4- a.c. -Stark shift of the ground state energy 

Since the efficiency of production of two-mode 
squeezed light is based on the resonant enhancement of 
two-photon processes, it is important to consider system- 
atically radiative corrections to the resonance frequencies 
in the implementation of the dynamics in Sec. [Til There- 
fore, we now evaluate corrections to the energy of state 
\tpi,m), Eq. (|27[). due to far-off resonance coupling in the 
limit of very small cavity decay rates. When considering 
the a.c. -Stark shift of state \ipi >m ), we find three contribu- 
tions, each associated to a different kind of coupling: (i) 
the a.c. -Stark shift due to the off-resonant laser coupling 
with the excited state at zero order in the mechanical ef- 
fects, Sujq ~ ft 2 / A for |A| ^> 7. It leads to a shift Suq of 
the dipole resonance frequency. The mechanical effects 
of the laser on the atoms give rise to (ii) a contribution 
which is linear in the number of vibrational excitation, 
and can hence be considered a renormalization of the trap 
frequency. This a.c. -Stark shift reads 



8v b w t] 2 cos 2 9 L il 2 tfb 
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2A(A 2 - v 2 + 7 2 /4) A 

(A 2 ~ v 2 + 7 2 /4) 2 + 7V 2 ~ A 2 + 7 2 /4 

Finally, off-resonant coupling of the cavity mode with 
the dipole transition gives rise to an a.c. -Stark shift of 
the cavity mode levels, which reads at leading order 
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(32) 



These shifts are in general not small and should be taken 
into account, when aiming at the resonant enhancement 
of certain processes over others. It should be remarked 
that the correction to Suij in Eq. (|3"2"]) which is at second 
order in rj arises from the mechanical effects of the in- 
teraction between resonator and center-of-mass motion. 
This term is non-linear, as it is a shift which depends on 
the number of vibrational excitation, but is a negligible 
contribution to Sujj. On the other hand, this term gives 
rise to an additional contribution to the a.c. -Stark shift 
of the center-of-mass motion, which is of the same order 
as Svb and depends on the number of photons. Its effect 
is detrimental, as the resulting spectrum of the center- 
of-mass excitations deviates from the one of a harmonic 
oscillator. In the system we consider we will neglect 
this contribution, focussing onto the regime in which the 
mechanical effects of the cavity mode can be neglected. 
This corresponds to situations, where the motion, for in- 
stance, is almost orthogonal to the cavity wave vector, 
I cos6U < 1. 



IV. SPECTRUM OF LIGHT AT THE CAVITY 
OUTPUT 

In this section we evaluate the spectrum of the light 
transmitted by the cavity mirror. The spectrum is best 
evaluated using the quantum Langevin equations for the 
operators Oj, at and b. The equations we obtain are 
rather involved, however the physical meaning of each 
term can be identified by comparison with the rates of 
the scattering processes discussed in the previous section. 



A. Quantum Langevin Equations 

We shall study the dynamics using the quantum 
Langevin equations (QLE) of the system. For conve- 
nience, we write the interaction Hamiltonian of the atom 
with the laser and the cavity fields as 



Hi, 



H al 



H n 



fii^B + aB^) , 



(33) 



where 



(31) 



B = n ( 1 - cos 2 6^(2^6+ 1) ) + irifL cos 6 L (tf +b) 
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3 = 1,2 \ 

—77 cos 9 c gj sin 4>j a j (b + &)j 
The QLE read 



cos" 



c (26 t 6+ 1) 



(34) 



ai(t) 
02 (*) 
6(t) 

<7(t) 

*«(*) 



i£iai(t) + icr(i) [-B(*) + , a x (t)] 

-Kioi(t) + v^afC*)) ( 35 ) 
2 a!(t) + «r(t) [B{t)\a 2 {t)} 

-K 2 a 2 (t) + V2^2~a l 2 n (t), (36) 
— *i/b(t) + ia(t) [B{t)\ b(t)] + icr(tf [B(t),b(t)} 
-K b b{t) + V^bb m {t), (37) 
~iA - 2] a(t) + (7,(t) (iB(t) + v^/ in (t)) , (38) 
2ia{t)B{tf - 2ia\t)B{t) - 7 [a-,(*) + 1] 
-2<rt(i)V7/ i "W - 2 ( x(i)V7/ m (0 t , (39) 



where cr z = er^cr — a<r\ and we have introduced the vac- 
uum input noises a™(i) (j = 1,2) of the cavity modes 
with corresponding decay rate Kj , the spontaneous emis- 
sion noise f m (t) at rate 7, and we also added a phe- 
nomenological input noise b m (t) acting on the atom's mo- 
tion, describing the heating at rate Kb due to the fluctu- 
ations of the trap potential. These four noise sources are 
mutually uncorrclatcd and have zero mean, while their 
second-order correlations have the form 

(a?(t)aT (f) f ) = (a 2 n (t)a?(t>)l) = S(t H), (40) 



(f in (t)f in {t^)=S(t-t% 
(b in (t)b in {t')i) = (N + l)S(t 
(b m {t)^b m (t')) =N5(t-t'), 



t'), 



(41) 
(42) 
(43) 



where TV is mean thermal vibrational number of the ef- 
fective thermal reservoir coupling to the atom center-of- 
mass motion [52j |. 

We assume that the laser is red-detuned and far-off res- 
onance from the atomic transition, i.e., A is negative and 
|A| is much larger than all the other parameters. This 
allows us to eliminate adiabatically the atomic internal 
degrees of freedom, and to assume that the atom always 
remains in the ground state \g), that is, a z (t) » — 1. 
Therefore we neglect the time evolution of a z , Eq. (|39[) . 
while Eq. ([38]) becomes 

a(t) = -(l~ iA) a(t) W(t) - Jlf n {t), (44) 

whose formal solution is 

a(t) = e~(^ A )V(0) (45) 
d se -(*- lA ) s [iB{t -s) + ^ff m (t - a)] . 



We now insert solution ([43]) into the other QLE and ne- 
glect the transient term because we are interested in the 
dynamics at times which are much larger than 1/|A|. We 
obtain 



Ol (t) 


= i5idi(t) + 


/ dse 
Jo 


(i-iA)s [B(t-s) - 


iVlf in (t-s)] [B(t)\ ai (t)} 


- n\ai(t) + v 




aa(t) 


= iS 2 ai(t) + 


/ dse - 
Jo 


(i-iA) s [ S ( t _ s )_ 


i^f n {t-s)] [B{t)\a 2 {t)} 


- n 2 a 2 (t) + v 




b(t) 


= -wb(t) + 


/ dse~ 
Jo 


(3-iA) s 


H/7/ m (*-*)] [B{t)\b{t)} 








— / dse" 
Jo 


-(i+iA)« 




(t-s)t] [B(t),b(t)]-K b b(t)- 


Ys/2h~ h b in {t), 





(46) 
(47) 

(48) 



r 



where we have not taken care of operator ordering, since, 
as we shall see, within the validity limit of our treatment 
these integral terms will generate only linear contribu- 
tions. 

At this point, we choose the laser angular frequency 
Wl so that 

61 =v' ; 5 2 = —v' 

namely, the laser frequency is tuned symmetrically be- 
tween the mode frequencies, which are spaced by a quan- 
tity 1v' . The angular frequency v' ~ v, and takes into 



account the a. c. -Stark shifts due to the mechanical cou- 
pling with laser and cavity modes, see Sec. lIII B~4l so that 
the two cavity modes are resonant with the motional side- 
bands of the laser light. Together with this choice of the 
laser frequency, we assume that the motional sidebands 
are well resolved, that is, v ^> |«fy|,0, Kj. 

In order to identify the resonant process, we move to 
a frame rotating at the effective vibrational angular fre- 
quency v' ~ v, (which has to be determined by solving 
the QLE) and we will neglect in the QLE all the terms 
oscillating at v' or larger. Denoting the slowly vary- 
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ing quantities by a\(t) = e ll/ *aj(i), a 2 (t) = e lv ' '* a 2 (t) , 
b(t) = e lv t b(t), after explicitly evaluating the commuta- 



tors we obtain 



a\(t) = i (1/ - <Ji) s{(t) + [ dse-^ +lA > \B{t- s^e 1 "'* +i^jf m {t- s^e™ 



a 2 ( 



- 9l cos 0i ( 1 - ^- cos 2 9 c (2b% + 1) j + r, gi sin 0i cos C (fe(i)e"^'* + 6 f (i)e w '*) 
(t) = i (i/ + fc) S2(t) + / dse-(*- 4A ) s \B(t - s)e lv ' 1 - i^f n (t - s)e lu ' 



-//■» c os02 ( 1 — 2" cos2 ^(2^6 + 1) ) + rjg:, mik >o cos f 



(5(t) ( 



6(t) = * (i/ - i/) 6(t) + / dse-(*-' iA ) s - s)e iv '' - i^yf n (t - s)e w ' 



Kia\(t) 



n 2 a 2 (t) 



i?7$!* cos6*l + 7?<7i sin 0i cos9 c a\(t)e %v * + 7752 sin 02 cos9 c a\(t)e 11 ' * 

dse^ +lA ) s - s)V Vt + i^/lf m {t - s^e™' 1 
— ir]Q cos 9 L + 7731 sin0i cos9 c ai(t)e w * + 77(72 sin0 2 cos 9 c a 2 (t)e~ w 1 — n b b(t) + y/2n b b m (t), 

I 



(49) 
/2^af(t)t, 

(50) 

/2^aj B (*)> 

(51) 

(52) 
(53) 



where we have introduced the noise operators a™(£) = making the Markovian approximation exp{— (7/2±iA + 

e -^'t a |n(t) ) a\ n (t) = e iy '*4 n (i), and 6 m (i) = e iv>t b in (t), imv')s} w (5(s)/( 7 /2 ± zA + imz/), for ™ = -1,0, 1. 

which are still delta-correlated. After long, but straightforward calculations we get the 

We insert in these equations the explicit expression final, effective QLE at leading order in the Lamb-Dicke 

for B(t - s), thereby neglecting the terms oscillating at parameter, which read 
v 1 or faster. We finally perform the time integrals by 



a 2 {t) 
b(t) 



i (1/ - <5i) a\(t) + xfb(t) - (ki + kil - iSi L ) a\(t) + V^afitf + V2R 1L ~^ L (t)^ + F u (54) 
i (v 1 + S 2 ) d 2 (t) + X 2~b(t) - (k 2 + k 2L + iS 2L ) 02(f) + V^a 2 n (t) + y/2~R 2L af L (t) + F 2 , (55) 
i (v' - v) b(t) + X16K*) - X 2 a 2 (t) - (Kb + ^2b - «i6 + ifo) &(*) (56) 
+V2k~bb m (t) + V2R 2b a 2 n L (t) - y/2K lb af L (tf + F b . 

I 



Let us now discuss each term appearing in the equations. 
The coupling coefficients are given by 



X2 = V^g'2 C0S <t>2 



COS Ul 



+ 



i tan cf> 2 cos ( 



A + t/-i7/2 A + i7/2 



Xi = rjftgl cos 0i 

X2 = V^92 C0S 02 

Xi = V^9i cos 0i 



cosOl i tan 0i cos ( 



A -1/' + 17/2 


A- 


hi 7 /2 


cos 9l 


itant 


b 2 cos 6 C 


A + v> + i 7 /2 


" A- 


Hi7/2 


cos 9l 


itan( 


b\ cos C 



(59) 

)• 

(60) 



(57) 

)• 

(58) 



A-v'-vy/2 A + i7/2 



and correspond to the Raman processes, in which laser 
photons are scattered into the cavity mode with a change 
in the center-of-mass excitation, see Sec. IIIIB 21 

New fluctuation-dissipation sources appear in the 
equations. We first discuss noise terms appearing in 
Eqs. (p)4"|) and In addition to cavity decay with 

rates Kj we find processes described by the decay terms 
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with rate k\l and K2L, and the corresponding Langevin 
noises a™ L (t) and 5 2 z,(*)> where 



7 |gi| 2 cos 2 0i 

KlL = 77 



K-2L 



and 



2 7 2 /4+(A-j/) 2 ' 
7 |g 2 | 2 cos 2 2 
2 7 2 /4+(A + z/) 2 ' 



—iv't 



< L (t) = f m (t)e 
(t) = f in (t)e il,,t 



-in 
a 2L 



(61) 
(62) 



(63) 
(64) 



These noises describe input-output processes between the 
cavity modes and external modes, mediated by the atom. 
They possess the same correlation functions of the spon- 
taneous emission noise f m (t), and at the timescales of 
interest, v't ^> 1, they are uncorrelated from each other, 
thanks to the oscillating factors. Note that 



KlL 



K2L 



7 gi cos ( 



2 7 /2 + i(A 



l 2. 92 COS( P2 



2 7 /2-i(A 



(65) 
(66) 



with KjL — \KjL\ 2 - They originate from the scatter- 
ing processes in which cavity photons are lost because 
they are absorbed and then spontaneously emitted by 
the atom, as has been discussed in Sec. MI B 31 

The noise and dissipation terms in Eq. (|56p . in addi- 
tion to the noise terms of the trap, are described by the 
decay terms with rate Kib and K2b, and the corresponding 
Langevin noise operators (i) and These pro- 

cesses originate from incoherent emission or absorption 
of a vibrational quantum accompanied by absorption and 
subsequent spontaneous emission of a laser photon. The 
emission of vibrational quanta takes place at rate 



K 2b = 77 



7 77 |ft| cos 9l 



2 7 2 /4+(A + i;') 2 



(67) 



while the rate of incoherent absorption of vibrational 
quanta is given by 



7 77 2 |ft| 2 cos 2 6> L 
2 7 2 /4+(A-i/) 2 ' 



(68) 



In particular, when A < then K 2 fc > K ib and the motion 
is cooled. Moreover, 



Klb = 



K2b 



7 rjfl cos8l 
2 7 /2 + i(A-z/)' 
7 77ft* cos Ol 
2 7 /2-i(A + j/)' 



(69) 
(70) 



with k 



\Kjb\ 2 . If we consider the dynamics described 
by these terms only, these incoherent phonon absorp- 
tion and emission processes lead to thermalization of the 



atomic motion at rate K26 — Kibj to a final effective mean 
vibrational number ruu = ^\b/{ K 2b — K ib) — |A|/4z/, as 
in standard cooling [48]. However, the noise associated 
with these incoherent phonon absorptions and emissions 
is correlated with the noise terms a™ L (t) and 5|" it) de- 
scribing scattering of cavity photons, because all these 
processes ultimately originate from spontaneous emis- 
sion. This is why the noise terms in the Langevin equa- 
tion for the atomic motion are directly expressed in terms 
of a % i L {t) and a 2 n L (t), making therefore this correlation 
evident. 

The operators Fj in Eqs. ([49|) -([53 |) represent non-linear 
terms, which describe the noise associated with the in- 
coherent part of the scattering processes discussed in 
Sec. MI B 21 These terms can be neglected with respect 
to the coherent processes, provided that 7 <C |A| and 
7 <C v. In particular, the second inequality ensures that 
rates xi and \2 differ appreciably, such that entangle- 
ment between the cavity modes can be established in a 
finite time [26]. We will focus on this regime, 7 -C v, in 
which we can thus neglect Fj in the effective QLE when 
evaluating the spectrum of squeezing. 

Finally, the frequency shifts of the two cavity modes 
and of the vibrational motion read 



SlL = 
(>2L = 



(A-E/)| gl | 2 COS 2 0! 

7 2 /4+(A-j/) 2 
_ (A + t/)|g 2 | 2 cos 2 <fe 

7 2 /4+ (A + i/) 2 
2Ar? 2 |r!| 2 cos 2 9 L ( 7 2 /4 + A 2 - v' 2 ) 



( 7 2 /4 + A 2 - v' 2 Y + z/ 2 7 2 
r? 2 1 ft 1 2 cos 2 <9 L - A 



(71) 
(72) 
(73) 



' A 2 + 7 2 /4 

and from their form one can recognize the a. c. -Stark 



shifts reported in Sec. MI B 4| with Sujj = 5 



Eq. (J32j>, and 8v b = 5 b tfb, Eq. (J3TJ), where now v — > v . 
Note that we have omitted a non-linear shift at second or- 
der in the Lamb-Dicke parameter, which affects both cav- 
ity modes and motion. As discussed in Sec. MI B~4l this is 
a small correction to 5jL, as it scales with rf , while it may 
have a relevant effect on the center-of-mass dynamics. It 
can be neglected in the limit ft cos 2 6l 3> 9j cos 2 8 C . Un- 
der this assumption, which we will consider in the rest 
of this manuscript, the spectrum of the center-of-mass is 
the spectrum of a harmonic oscillator, characterized by 
equidistant energy levels. 

As the dynamics we seek relies on resonant interaction 
between the cavity modes and the vibrational motion, the 
two cavity modes should be exactly at resonance with the 
sidebands of the driving laser. Equation (|73[) provides 
an implicit equation for the actual vibrational angular 
frequency v' . In the parameter regime ?7|ft| <C |A|, v we 
find with good approximation 

2A?7 2 |ft| 2 cos 2 9 L ( 7 2 /4 + A 2 - v 2 ) 



( 7 2 /4- 



v 2 ) + z/ 2 7 2 
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-^IW ^^ . (74) 

Taking also into account the frequency shifts of Eqs. (jTTj) - 
(|72p. the resonance conditions are finally 

= S 1L + v' (75) 
S 2 = S 2L - v'. (76) 



In the parameter regime 7 <C using conditions (|75|) - 
([76]) . we arrive therefore to the final QLE, describing the 
coherent interaction between the two cavity modes and 
the vibrational motion, competing with losses and noise 
processes due to spontaneous emission, cavity decay, and 
vibrational heating, 



a\(t) = xW) - (*i + Mi)aJ(i) + v^Ta^tjt + V2R 1L a™ L (tf , (77) 
02(f) = Xa&(*) - («a + ^2i) oa(t) + \/2^af(i) + V2R 2L a^ n L {t) , (78) 
6(t) = Xia\{t) - x* 2 a 2 {t) - (/s 6 + K 2fc - re 16 ) 6(t) + V2k~ b b m (t) + V2R 2b a% l L (t) - v^ftiiO&Ct)* . (79) 



B. Evaluation of the spectrum of squeezing 

We now use Eqs. l[77|) - ([79"j) in order to determine the 
stationary spectrum of squeezing of the light at the cavity 
output. We consider 



I mt (t) = al ut {t) +a°"*(t) t -al u \t) -a° 2 ut {t)\ (80) 
Il ut (t) = -ilal^i^-al^^ +a^ ut {t)-a^ ut (t)\}l) 



corresponding respectively to the difference between the 
amplitude quadratures, and the sum of the phase quadra- 
tures of the two sideband modes. These are the quadra- 
tures exhibiting two-mode squeezing in the case of pulsed 
excitation in this setup, see [25l . [2q |. The output cav- 
ity fields a° ut (t) in Eqs. ([5D ]) -(|5T ]) are given by the usual 
input-output relation 



'(t) = y/2k~ aj (t) - af(t), j = l,2. 



(82) 



The spectrum of squeezing can be calculated by evaluat- 
ing the Fourier transforms 



dte ^tjout 



(t), 



(83) 



and using the fact that at the stationary state it is 



(I. 

where we have normalized the spectrum so that the shot 
noise level corresponds to S±(lo) = 1. Two-mode squeez- 
ing is found when one spectrum of squeezing takes val- 
ues below the shot noise limit at some u. From the 
Fourier transform of Eqs. (|77jl - l|79jl one can see that 
S+(ui) = S—(u>) = S(lo), which implies that in the 
present case two-mode squeezing is equivalent to EPR- 
like entanglement between the two output cavity modes. 
This is easily verified by applying a sufficient criterion for 
entanglement, such as the "sum" criterion of Duan et al. 
f53[ | , or the product criterion of Ref . [H, HE] • With the 



chosen normalization for the output cavity modes at u), 
the sum criterion reads 



S+H + S_H < 2, 
while the product criterion gives 

< i, 



(85) 



(86) 



so that in our case both criteria imply that the two out- 
put modes are EPR-likc entangled as soon as S(oj) < 1. 
The squeezing spectrum S(u>) can be obtained from the 
Fourier transform of the Langevin equations after long 
but straightforward algebra, yielding a cumbersome ex- 
pression which will not be reported here. This expression 
becomes considerably simpler in the limit |0|, \gj\, 7 <g; 
|A| and r] <C 1. In this limit the additional loss pro- 
cesses due to spontaneous emission, associated with the 
rates Kji, and Kj b (j — 1,2), are typically negligible, that 
is, KjL, Kjb <C k. Moreover, we consider the case of ion 
traps, where heating of the atomic motion is negligible 
with respect to all radiative noise sources (56|. Finally, 
as the two cavity modes are very close in frequency, they 
will have very similar properties, in particular we can 
take k± — k 2 = k. In this parameter regime the main 
aspects of the squeezing spectrum can be grasped from 
its analytical expression. One finds 



S(oj) = 1 - 



k 2 (e 4 - e 4 ) 



(K 2 +UJ 2 ) {U 2 -9 2 ) 2 +UJ 2 K 2 



(87) 



where = ^J\x 2 \ 2 — |xi| 2 as given in Eq. and 



IX2| 2 + |xi| 2 



2xiX2 



and we have used that Xj — Xj when 7 <C |A| (see 
Eqs. (|57)l - (|6^|) '). Note that due to the transformations 
which we have applied, the results which appear around 
uj = in S(u>) describe quantum correlations of noise 
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components in the optical signal at lul — is' ± to with 
those at u L + v' ± w, i.e. correlated fluctuations of the 
two modes at the same offset from their center frequen- 
cies. 

From Eq. (|8T[) one notes that the properties of the 
spectrum are mainly determined by the ratio Q/k. In 
fact, the denominator in Eq. ([57)) has always three 
poles in the lower complex half-plane, one which is al- 
ways imaginary at loq — —in, and two poles at tu± = 
-in/2 ± ^/0 2 - k 2 /4. Therefore, when Q/k > 1 the 
two poles at uj± have a nonzero real part and the spec- 
trum is characterized by three well-separated inverted 
Lorentzian peaks, one at u> = with width (FWHM) 
2k and the other two symmetrically placed atw» ±8, 
with FWHM k. At the center of these peaks one has 
S = (S/6) 4 ~ [(|x 2 | 7 IXi|)/(IX2| + |xi|)] 2 , approach- 
ing S(u>) = 0, i.e., infinite two-mode squeezing, for 
\X2\ — \xi\- Therefore, when > k we find two-mode 
squeezing within three narrow bandwidths around u> = 
and u> = ±0. In the opposite case of < k/2, the 
three poles are all on the imaginary axis, and the spec- 
trum shows only one inverted Lorentzian peak at u> = 0. 
When Q/k -C 1, this peak becomes very narrow, with 
FWHM - 2& 2 /k. 

It is remarkable that even for k > one finds almost 
perfect squeezing in the difference of amplitude quadra- 
tures at lj — 0. This can be understood, as in the regime 
we consider the scattered photons due to spontaneous 
emission are negligible with respect to those lost through 
the output cavity mirror (KjL,Kjb <C k). This implies 
that most of the intracavity photons are detected at the 
output. These photons are almost perfectly correlated at 
ijj = and therefore would give S(0) ~ 0. In this regime, 
a large cavity decay rate k has only the effect of narrowing 
the squeezing bandwidth. On the contrary, if the photon 
scattering by spontaneous emission is no more negligible, 
two-mode s quee zing soon degrades, even at u> = (see 
for example [32|). 

The presence of the three-pole structure in the squeez- 
ing spectrum is novel with respect to the spectral features 
usually encountered in the parametric oscillator (either 
below and above threshold, see e.g. Ref. [Hj]). This struc- 
ture is due to the coherent interaction of the two cavity 
modes with the quantized atomic motion, i.e. it arises 
from the coherent microscopic processes underlying the 
dynamic establishment of quantum correlations. The pe- 
culiar spectral properties can be exploited to achieve opti- 
mal broadband two-mode squeezing when the three peaks 
merge, which happens when = k. In this case one 
easily sees from Eq. (|5Tj) that 



5(w) = 1 - 



> ( Q 4 _ 

Q e +w 6 



(89) 



i.e., one has large, uniform squeezing for a wide band- 
width of frequencies. The fact that = k is the condition 
for the best two-mode squeezing in the output can be eas- 
ily understood noticing that is the angular frequency 
at which two-mode squeezing inside the cavity is periodi- 



cally built up (see Ref. [25|, [26|). Therefore, when = k, 
squeezing is generated inside the cavity at the same rate 
at which it is transferred to the output field. In contrast, 
in the other two cases, squeezing is not efficiently gen- 
erated in the output because either the output coupling 
happens too fast, i.e. before full intra-cavity squeezing 
is established, or the output coupling is sufficiently slow 
to allow that energy is stored inside the cavity and the 
two-mode squeezing is coherently re-converted into inde- 
pendent states before it is coupled out. 



C. Results 

We now consider the exact squeezing spectrum S(lu), 
defined in Eq. ([84)1 . where the operator ([80)1 is evaluated 
using the output relation ([B"2")l , such that the operators 
Oj(i), a](t) are the solutions of the QLE ([TTJ-JTS]) includ- 
ing all noise and loss terms. 

The parameter regime we consider has been discussed 
in detail in [2(| . We take a F= « F = l atomic tran- 
sition with the quantization axis B along the cavity axis, 
and B, kL, and laser polarization El mutually orthogo- 
nal. Relation 7 <C v can be fulfilled by the intercombina- 
tion line of an Indium ion at 7 = 2~k x 360 kHz in a trap 
at v = 2tt x 3 MHz, for which r\ ~ 0.1. We consider 
a geometrical configuration corresponding to 0l = 0, 
9 C = 7r/2, and 4>i — <p2 = 0, which means that the ion 
motion takes place along the direction of the laser beam 
and orthogonally to the cavity axis, and that the trap 
center coincides with an antinode of the cavity modes. In 
such a case, the coupling coefficients of Eqs. ([57l) - ([60l) are 
determined by the first term only. Moreover in this case 
Dcos 2 6>l ^> gj cos 2 C ~ 0, and therefore, as discussed 
in Sees. IHI B 41 and IIV Al the small Stark-shift correc- 
tion to SjL, scaling with if can be neglected. If we con- 
sider that the ion couples to two non-degenerate polariza- 
tion modes of a resonator with vacuum Rabi couplings 
g ps 2ir x 0.6 MHz, and we take laser Rabi frequency 
tt = 2tt x 18 MHz and detuning A = 2tt x 60 MHz, we 
obtain Q/2n ~ 7.9 kHz, see [26|. The condition 3> k 
is found for a finesse T ~ 10 6 and free spectral range 
5lj = 2n x 1 GHz, so that k — 2n x 1 KHz. The cor- 
responding spectrum of squeezing is displayed in Fig. [3] 
(full line). It exhibits three minima at u> — 0, ±0, which 
correspond to three separated regions of narrow-band 
squeezing. The two bands around u> = ±0 have width 
k, while the central one has width 2k and shows almost 
100% squeezing. These features are well reproduced by 
the analytical expression ([57)1 (see dashed line in Fig. [3]), 
except that the latter predicts very large squeezing also 
for the peaks at uj = ±0. The success of the simplified 
expression ([57)1 is due to the fact that, with the cho- 
sen parameter values, the loss rates due to the various 
scattering processes are at least ten times smaller than 
the cavity output loss rates k\ — «2 = k, and there- 
fore do not have a relevant effect on the spectrum. We 
have also considered a realistic ion vibrational heating 



rate Kh = KbN — 2ir x 0.1 kHz, which however gives an 
effect which is negligible even with respect to that due to 
photon scattering. The appearance of three minima is a 
novel behavior to our knowledge, and it arises from the 
coherent microscopic dynamics, as modulates the ex- 
change of excitations and correlations between the cavity 
modes and the center-of-mass motion. 




FIG. 3: Squeezing spectrum S(to) as a function of the side- 
band frequency u> (in units of O) when Q ~ 8k. Its behavior 
is well reproduced by the approximate analytical expression 
of Eq. (|87p (dashed line). The parameters are discussed in 
the text. 



The most interesting regime of broadband two-mode 
squeezing, when ~ k, is shown in Fig.[4l which displays 
the squeezing spectrum in the case of the same parameter 
values of Fig. [3] except for a lower cavity finesse, T ~ 10 5 , 
implying n = 2ir x 10 kHz. The three minima merge into 
a single broad one, centered around u> = 0, whose width is 
determined by n = 0. Also in this case one gets almost 
perfect squeezing at the center, and these features are 
well reproduced by the simple analytical expression of 
Eq. fg3]) (dashed line in Fig.gJ). 
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FIG. 4: Squeezing spectrum S(u) when Q = k. Its behavior 
is well reproduced by the approximate analytical expression 
of Eq. (|87|) (dashed line). The parameters are discussed in 
the text. 




FIG. 5: Squeezing spectrum S(u) when ~ k/10. Its behav- 
ior is well reproduced by the approximate analytical expres- 
sion of Eq. (|87[1 (dashed line). The parameters are discussed 
in the text. 



V. CONCLUSIONS 



Finally, in Fig. Owe consider the case k > 0. We have 
still kept the parameter values of Fig. [3l but we have 
now considered a cavity with finesse T ~ 10 4 , imply- 
ing k == 2tt x 100 kHz. The squeezing features are visi- 
bly worsened, as in this regime losses are faster than the 
typical time scale in which correlations between the field 
modes are established. One has still two-mode squeezing 
around lu — 0, but with a very narrow bandwidth which 
is roughly given by Q 2 /2k. 



The resonance fluorescence of a confined, single, laser- 
driven atom exhibits EPR-entanglement, or two- mode 
squeezing, in the field modes which interact resonantly 
with the Stokes and anti-Stokes transitions created by 
the atomic motion. By coupling these sidebands to a 
high-finesse optical cavity, we have shown how to cre- 
ate continuous-wave (cw) two-mode squeezed light out- 
put from the cavity. At the microscopic level, the pro- 
cess is based on the mechanical effect of light, which 
allows for quantum-coherent generation and control of 
entanglement between the motion and the cavity modes. 
The scattering processes have been characterized and de- 
scribed in simple physical pictures using scattering ma- 
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trix theory, and the squeezing spectrum has been calcu- 
lated using Quantum Langevin Equations. 

Peculiar novel spectral properties are predicted for the 
squeezing spectrum of the cavity output. They may be 
divided into three regimes of the cavity output rate k rel- 
ative to the frequency of creation of two-mode squeez- 
ing inside the cavity. The squeezing spectrum can consist 
of a single peak (n > 8), three peaks (k < ©), or one 
broad, homogeneous band (k = O). Simple analytical 
approximations have been derived for the three relevant 
regimes. 

The squeezing spectrum in the different parameter 
regimes has been calculated for an experimentally ac- 
cessible case of a single trapped ion as a specific ex- 
ample. The results for this specific system show all 
the features predicted by the general derivations, ex- 
hibiting novel spectral properties of two-mode squeez- 
ing which are novel when compared with conventional 
Optical-parameter-amplifier-type of sources. 

In particular, starting from the most fundamental in- 
dividual quantum systems, a single atom and an optical 
cavity, we have designed a nonlinear optical source. This 



is therefore a paradigmatic model system exhibiting the 
connection between microscopic, quantum-coherent dy- 
namics and macroscopic nonlinear device properties. Its 
efficiency and the high-degree of control one can achieve 
on its dynamics offer promising perspectives for the re- 
alization of quantum light sources for quantum network- 
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